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r D ■ Abstract 

Correlation functions and form factors in vertex models or spin chains are known to satisfy 
certain difference equations called the quantum Knizhnik-Zamolodchikov equations. We find 
similar difference equations for the case of semi-infinite spin chain systems with integrable 
. boundary conditions. We derive these equations using the properties of the vertex operators 

and the boundary vacuum state, or alternatively through corner transfer matrix arguments for 
the 8-vertex model with a boundary. The spontaneous boundary magnetization is found by 
solving such difference equations. The boundary S-matrix is also proposed and compared, in 

o 
in 



X 



the sine-Gordon limit, with Ghoshal-Zamolodchikov's result. The axioms satisfied by the form 



Q\ , factors in the boundary theory are formulated. 



1 Introduction 



Following the success of conformal field theory in two dimensions, much attention has recently 
been focused on trying to understand the algebraic structure of solvable lattice models and massive 
integrable quantum field theories |j],||.|3],|j . Certain quantum deformations of affine Lie algebras play 
a central role in the description of these systems. One of the remarkable results of these studies is 
the discovery of the difference analogues of the Knizhnik-Zamolodchikov (KZ) equation ||. These 
'quantum KZ equations' are satisfied by both the correlation functions and the form factors of 
integrable models |I],|2|,Q||. They allow us to study the off-shell properties of the models. 

In this paper we establish similar results for the XYZ Hamiltonian on a semi-infinite chain, with 
an interaction at the boundary corresponding to a magnetic field in the z-direction: 

-. oo 

R b = ~ 9 E ((! + r K+K + (i - rK+i^ + Aoi+io*) + hoi (i-i) 
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Integrable models with special boundary conditions which preserve integrability, have been 
studied both in lattice and continuum theories, e.g., (7]||. In general frameworks, the boundary 
interaction is specified by a reflection matrix K for lattice systems M, or by a boundary 5-matrix 



for quantum field theories [10|. Integrability is guaranteed by the fact that the they satisfy the 
boundary Yang-Baxter equation. The Hamiltonian ( |1.1| ) is obtained from a commuting transfer 
matrix constructed by such a iT-matrix This model is particularly interesting because the 

sine-Gordon model |1(J is its continuum limit. 

In a recent paper [O], it is shown that the space of states of the XXZ spin chain with a 
boundary can be described in terms of the vertex operators $ e (C) 5 **(C) associated with the "bulk" 
(infinite chain) XXZ model j2|. Explicit expressions for the two vacuum states \i)bi {i = 0,1) of 
the boundary XXZ Hamiltonian are obtained by using the bosonization formula for the vertex 
operators, and then they are used to obtain the boundary excited energy and the boundary S- 
matrix. 

It is an interesting problem to extend the analysis of the off-shell quantities such as the corre- 
lation functions and the form factors in the bulk theories to the boundary cases. In this paper, we 



show that the boundary analogues [14] of the quantum KZ equation are satisfied by these quanti- 
ties in the XXZ and XYZ models. For the XXZ model, our results follow from the results in fl3| . 
Assuming similar properties of the vertex operators and the vacuum states, we can extend the 
results to the XYZ model. Alternatively, for the correlation functions we can use corner transfer 
matrix arguments |16|,]17| to derive the same difference equation. 

In Section 2, we will derive the ^-difference equations satisfied by the correlation functions of 
the XXZ and XYZ spin chains with a boundary. In Section 3, we derive the same equations for 
the XYZ chain using a graphical corner transfer matrix argument similar to that used in [l7|| . In 
Section 4 we discuss the boundary magnetization (the vacuum expectation value of the boundary 
spin operator erf) of the XYZ model. In the case where h = 0, one can solve the corresponding 
difference equation, and the result is simply the square of the bulk magnetization of the eight- vertex 
model. A similar phenomenon was observed in the case of the XXZ model with a boundary J13| . 
We conjecture the form of the boundary S-matrix in Section 5. Lacking bosonization formulas for 
the vertex operators of the XYZ model, we cannot derive it explicitly. However, our result reduces 
to the appropriate sine-Gordon boundary S'-matrix of []l(| in the scaling limit. Finally, in Section 
6, we present the analog of Smirnov's axioms § for the form factors. 



2 The boundary (/-difference equation 

In this section, we derive the (/-difference equation satisfied by the boundary n-point function, i.e. 
the expectation value of the product of type I vertex operators between the vacuum states of the 
boundary Hamiltonian. This equation was first discovered by Cherednik [Q. We will derive this 
equation first for the boundary XXZ model in the formulation of jljj, and then extend it to the 
boundary XYZ model. 
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2.1 The XXZ model 



The boundary XXZ model was formulated in [13] in terms of the vertex operators of Uqi^slo). Let 



us recall the basic relations for type I vertex operators <J? £ (C) and 3>g(C) = &-e{— Q~ 1 (), £ = ± ||j 
and the boundary ground state \0) B of the XXZ Hamiltonian defined in | |13| |. For the definition of 
the i?-matrix R(C) and the boundary reflection matrix K(Q) we refer the reader to |TI| (see (A.l) 
and (A. 2), and (2.2) and (2.3) therein, and see also (A. 11) for the scalar g which appears below). 
The type I vertex operators commute as 

E «S;S(Ci/c 2 )^ 1 (Ci)^(c 2 ) = <mc 2 )<mci) , (2.1) 

f' f' 

and have the inversion property 

J*«(0*«(0 = W. (2-2) 

The boundary ground state and its dual have the following reflection properties with respect to 
type I vertex operators: 

E^(0*«'(C)io> B = * e (r 1 )|o>fl , (2.3) 

s' 

X)fl<oi^'(r 1 W(c) = fl(oi**(o • (2.4) 

e' 

The relation between these and the spin chain Hamiltonian will be discussed in the context of 
the XYZ model in Section 2.3. Here we note only that the (renormalized) XXZ Hamiltonian is 
given by H = const. dT B (() / 'd( ^, where the transfer matrix is defined by 

T B (C) =gY.K(C 1 )K!'(0<S>e'(C). (2-5) 

The ground state \0} B satisfies T B (()\0) B = \0) B . 
Consider the boundary n-point function, 

G((l,- •-,Cn)= E B(0|* E l(Cl)---* E n(Cn)|0)^ £l ®---(8^, (2.6) 

where u ± are the weight vectors in the 2-dimensional vector space V = Cv + Cv~ . Note that 
R(C) G End c (V(g)V) and K{Q) € End c (y). Equations (O), (E3) and (p|) imply 



Here, 



Rjj+l\Cj/ Cj+i) G(Ci, • • • , Cj, O+i) " ■ ■ j Cn) — -fj'i+iC(Ci) ' • • i 0+1) 0' ' ' ' ' Cn), 

^n(Cn) G(C 1 ,---,Cn)=G(Cl,---,C 1 ), 
^l(Cl) G(d _1 , • • • , Cn) = G^Cl, C2, • • • , Cn). 



P jj+ i(v £l ® • • • ® ® f £j ' ® • • • ® t; £n ) = t> £l ® ■ ■ ■ ® u £j ® v £j+1 ® • • • ® v £n , 
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R jj+1 (()(v £l ® ■■■ ®v £ ") = V R £ j' e J +1 (0(v £1 <8 ■■■®v £ 'i ®v e 'i+i 
Kj(()(v £l • • • (8) w £ ") = ^^/(C)(u £l (8) • • • <g> v e 'i (8 ■ ■ ■ <8 w £ "), 
Kj(()(v £l ®---®v £n ) = ^2K_ £ J .(-q~ 1 ()(v £l (8>---<8u £ J ®---<8u e "). 



(8 w £n ' 



Prom these equations we obtain 

x Rij(CiCj) ■ • ■ Rj-ij(.Cj-iCj)Rj+i ■ • • Rnj((n(j) x 

x K j (( j )R jn (( j /( n )---R jj+1 (( j /( j+1 )G(( 1 ,---,(n)- (2.7) 



This is a version of Cherednik's equation |14 |. 



2.2 Consistency conditions 



We have shown that the boundary n-point function • • • , £ n ) satisfies (2.7) if relations (2.1), 



and (2.4) are satisfied, without using any further details of the XXZ model. Therefore, if we 



have more general settings for these relations, we get more general solutions to equation fl2.7|) . We 
remark that these three relations along with (|2.2| ) imply the following consistency conditions for 
R{() and K(C): 

(i) Yang-Baxter equation 

i?12(Cl/C2)#13(Cl/C3)i?23(C 2 /C3) = ^(Ca/Cs)^ IstCl/Cs^Cl^)- (2-8) 

(ii) Unitarity relation 

i*i2(Ci/C 2 )i*2i (C2/C1) = id. (2.9) 

(hi) Crossing relation 

(iv) Boundary Yang-Baxter equation 

^2(C 2 ) J R2l(ClC2)^l(Cl)^12(Cl/C2) = i?2l(Cl/C2)^l(Cl) J Rl2(ClC2)^2(C2). (2.11) 

(v) Boundary unitarity relation 

KiOKiC 1 ) = id. (2.12) 

(vi) Boundary crossing relation 

<(0 = E flMWK'r 1 )- (2.13) 

a',b' 

For completeness, we add a few remarks following the results in The properties of the vertex 
operators entails the following relations for fl2.5|) : [7b(£i),Tb(C2)] = 0, Tg(l) = 1, Tb(C)?b(C X ) = 



4 



1, T B (-g- 1 C _1 ) = T B ((). Therefore, i£t(() is an eigenvalue of T B (C), wehavet(l) = 1, t(C)*(C _1 ) = 
1, and ti-q^C 1 ) = t((). A vector \v) (or satisfies £ e , (Q$ e >(£)\v) = i(C)^e(C _1 ) l^> (or 
2 e (v[$e(C~ 1 )-K'f / (0 = t{()( v \®* e >(0) if an d only if it is an eigenvector of T B (() with the eigenvalue 
i(£). In this situation, the n-point function 



G(Cl,-'-,Cn)= E (v\*eAtl)---*e n (Cn)\v)v £ * 



(8) ti £n 



(2.14) 



satisfies the same equation ( |2.7| ) because the effects of t(C) in if (0 and if (0 cancel. In particular, 
B<l|* ei (Cl)-"* e n(Cn)|l)B satisfies (|3) (see 3.1 in @). 



2.3 The XYZ model 

In this section we consider ( |2.7| ) in the context of the XYZ model, extending the setting of Section 
2.1. Unlike the XXZ case, we do not have a complete solution for the XYZ model, because we do 
not have the corresponding mathematical machinery, in particular, a bosonization scheme. Certain 
results can, however, be generalized. Following |]l5| we assume the existence of vertex operators. 
We further assume the existence of the boundary vacuum states, and derive ([O]) for the boundary 
n-point function. See, however, Section 3, where we give a physical argument that supports the 
existence of such a setting. 

To define the model, we use the variables 



P 



e k , 



7TA 

e ax . 



c 2 , 



(2.15) 



where K, K' are /, I' in Baxter's notation fil| , and A and u are as in (10.4.21) of 16]. The variable 
p is the elliptic nome, and p = corresponds to the XXZ case. We restrict our discussion to the 
principal regime, < p 1 ! 2 < — q < < 1- Let 



snh(ii) 
(z;qi, ■ 



isn(in), cnh(n) = cn(iu 
in 



dnh(n) = dn(iti), 



,q m )oo= I| (1 

ni,...,n m =0 



qr---q^ m z), Q q {z) 



(z;q) 00 (qz 1 ;q) 00 {q;q)c 



The elliptic functions sn(n) etc., are found in the appendix of [16]. Unless otherwise stated, the 
elliptic nome for these functions is p, and the corresponding modulus and conjugate modulus are 
denoted by k and k' . In ( 2.2SQ etc., we use elliptic functions with elliptic nome 



-ILL. 

e i 



(2.16) 



modulus ki and conjugate modulus k'j. 

We first give solutions R and K to the consistency conditions (i)-(vi). The solution R{Q) to 
equations (i)-(iii) for the XYZ model is 

MO d(C)\ 

6(0 c(0 

c(0 HO 

\d(0 a(C)J 



R(0 



(2.17) 
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where the Boltzmann weights are 

1 snh(A — u) , . 1 snh(u) 



«(0 



with 



MC) snh(A) 



c(C) 



M(C)snh(A)' v ^ MC) 

1 (p 2 ;p 2 )^e p2 (g 2 )9 p2 (^) 



snh(A — u)snh(u) 
MO snh(A) ! 



/"(C) k(C 2 ) (p;p); 



1 (q 4 z 1 ;p,q 4 )oo(q 2 z;p,q 4 ) 



Q p 2(q 2 z) 

>(pz~ l ;p, q 4 )oo(pq 2 z;p, <? 4 ) c 



(2 18) 

«(z) (q ,4 z;p,9 4 )oo(^ 2 ^ _1 ;p,9 4 )oo(^;p,9 4 )oo(w 2 ^ _1 ;p,q 4 )oo' 

Usually, will not denote the p and q dependence explicitly in cases when no confusion should arise 
as a result of their absence. 

The general 3-parameter family of solutions K(() to ( 2.1 1| ) was given in [11| and \12\. In the 
following, we restrict to the diagonal solution K(Q) with a single parameter 



7T 77 



Let 



where 



f(C,r) 



K((;r), K(C;r) 



snh(r]+u) 
snh(»?--u) 



(2.19) 
(2.20) 



/(C;r) 

ip(z;r) 



-1. 



<p(z;r) 
f(z 

(prz;p 2 ,q i ) O0 (p 2 r~ 1 z;p 2 ,q 4: ) O0 (rq i z;p 2 ,q 4: ) OD (pr~ 1 q 4: z;p 2 ,q i ) O0 
(prq 2 z;p 2 , q i ) OQ (p 2 r~ l q 2 z; p 2 , q A ) OQ (rq 2 z; p 2 , q 4 )^ (pr~ 1 q 2 z; p 2 ,q A )oo 



{q b z 2 ;p 2 , q*)oo{pq 2 z 2 ;p 2 , q s ) 



>z 2 ;p 2 ,q 8 ) OQ (p 2 q 2 z 2 ;p 2 ,q 8 ) 



(q 8 z 2 ;p 2 ,q 8 ) 00 (pq' i z 2 ;p 2 ,q 8 ) 2 X) (p 2 z 2 ;p 2 ,q 8 ) 00 

The scalar factor f((;r) is chosen so that (|2.12j ), ( |2.13| ) hold, and KZ{C) is analytic in the region 
—q < \r\ 1 ^ 2 < In terms of the additive parameters defined in ( 2.151 ) an d ( |2.19| ), the region we 
consider is given by < u < —rj<X<K'. 

Now let us recall the vertex operators defined for the elliptic algebra corresponding to the 
i?-matrix ( 2.17| ) in |l5|][]. The type I vertex operator 



acts on graded vector spaces 



with character 



$ e (C) : n {i) -» Ti^ 
n^ = ®f = ,nf (* = o,i) 

£dim7^V- 



d=0 



n^ = o(i-^ +1 )' 



1 In JlT]], $ e (C) was introduced satisfying (2.1), and (2.2) with a different normalization, i.e., g = 1. In this paper 
we follow the normalization given in J13]. 
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The vertex operator is expanded in the form 

where : H.y — > T^-n ■ ^ e normalizations are such that $-^(0) = [1), and $ +j o|l) = |0), 
where TLq = C\i). In this normalization, we find that Eq. ( |2,1| ) holds with the elliptic 72-matrix 
(gl7D , and Eq. Q holds with || 

= (pq 2 ; p) oo {p<f ; p, g 4 ) oo {pq 2 ; p, g 4 ) oo (g 2 ; g 4 ) oo , 2 21 s 

5 fep)oo(M 4 ;p,g 4 ) 2 o(g 4 ;g 4 )oo " 1 ' j 

We now turn to the boundary XYZ hamiltonian. In analogy with the XXZ case, we now assume 
that there exists a vacuum vector [0}b (resp. b(0|) in an appropriate completion of ©d : even?^ \ 
(resp. ©d:even^d° )' which satisfies fl2.3|) (resp. (|2.4j) ). In Section 3, we give a graphical argument 
for the construction of these states. However, since we do not have a bosonization formula for the 



vertex operators, we cannot give an explicit bosonic construction of this state, as we did in [13| for 
the XXZ Hamiltonian. 

In the XXZ model, we found that there exists another vacuum state, the boundary bound state 
|1)b, with an excitation energy 

A(C . r) mill - 1 ®ArC 2 )&A<l 2 rC- 2 ) 

ie|(C;r-i) ee q 4rC- 2 )&A<i 2r C 2 y { ' } 



It turns out that this formula remains valid without any change for the elliptic K (£) given by ( 2.20 ). 
Therefore, for \r\ > \q\, we assume that there exists a vector \1)b (resp. #(1|) in an appropriate 
completion of (B&evenMd : j (resp. (Bd-.even'H^*)' which satisfies 

^K*(C)*AQ\l) B =Afcr)* E (.C 1 )\l)B , (2-23) 

s' 

J2B(m>(C 1 )K!,(0=A(( ] r) B (l\m) ■ (2-24) 

e' 

Under these assumptions, we can formulate the boundary XYZ model in an analogous way to 
the boundary XXZ model. In particular, the correlation function 

E B{i\^e 1 (Cl)---^e n (Cn)\i)BV ei <g> ■ ■ ■ ® V £ " 
ei, ■■■,£„ 

with elliptic vertex operators satisfies the g-difference equation 

We define the transfer matrix of the boundary theory as in (|2JJ) with the appropriate elliptic 
generalizations of g, K and 3> e (C)- The renormalized form of the Hamiltonian Ql.ip is defined from 
this transfer matrix by 

rrrenor _ 7TSnh(A, fc) d 

4^ C d( TB{C) C =i- (2 - 25) 

We find that the parameters in equation Ql.l| ) ar e related to the elliptic parameters above by 

r = fcsnh 2 (A,£0, A = -cnh(A,A:)dnh(A,A;), h = -snh(A, k) Cnh(?7, fc) ■ (2.26) 
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This Hamiltonian is normalized so that its lowest eigenvalue is 0. From ( 2.22| ) and ( 2.25| ), the 
eigenvalues e^(r) (i = 0, 1) of Hb on are 



=(0) 



(r) 
(r) 



e(l)sn(-2I'(77 + iK)/A,^) 
e(l) 



if —1 < r < — |g| 
if |g| < r < 1. 



(2.27) 
(2.28) 



-ze 



„i0 



tysn(-2JVA, fcj) 
where we have defined the following function 

e(0 = -Lnh(A,A;)dn(— 0,*,), 

A 7T 

6(1) = -^snh(A, k). 

The reader should be careful about the nome of the elliptic functions. (See ( |2.16| ).) The function 
e(£) differs from the one in the XXZ chain 13] due to the difference of the normalization factor in 
front of the derivative of Tb{£) in ( 2.25| ). It reduces to the XXZ excitation energy when p = 0. 



The excited states are created by the action of type II vertex operators = =t) Ijj on 



the vacuum vectors \i)b- 



(2.29) 



The energy spectrum of the one particle excitation is evaluated as e'*'(r) + e(£). The excitation is 
therefore massive with mass e(l) fig]. 



3 The Corner Transfer Matrix 

In this section, we show how within the context of an inhomogeneous boundary vertex model it is 



possible to construct lattice realizations of < I ) e (C)) ^(C)> and #(0| that obey relations ( 2.1 
( |2.3f) and ( |2.4D . We then present a rather physical picture of the origin of the difference equations 
for correlation functions of local operators in this vertex model. 



3.1 The Partition Function 

Following Sklyanin Q, we build a lattice with 2M vertical and 2N horizontal lines as in Figure 
1. Here, in the vertical direction we impose the cyclic boundary condition. Let us call this lattice 
Cmn- Let the horizontal lines carry alternating spectral parameters £ and l/£, and the vertical 
lines £. We consider the region < p 1 / 2 < —q < \r\ 1 / 2 < |C/£| -1 , IC£I _1 < 1- We use the word 
ground state to refer to the lowest energy states in the /i = 0,r = 0,A— >— oo limit. There are two 
such ground states, one of which is shown in Figure 2. Note that if h ^ then the degeneracy of 
the ground state energies is resolved, and only one of the ground states corresponds to the lowest 
energy. Nevertheless, we call both of them ground states for all values of h. Corresponding to these 
two ground states, we have two antiferromagnetic sectors. Let us label these sectors by i = 0, 1. 
When we consider the partition function and correlation functions, we choose one of the ground 
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L/C 

c 

i/C 



nc 



2iV lines 



e e e 

2M lines 



Figure 1: The inhomogeneous 2M x 2iV lattice Cmn- The circles on the right hand side denote 
reflection matrices K(Q, and those on the left hand side denote K{— (? _1 C _1 )*- 
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Figure 3: The inhomogeneous lattice split into four sections 



state sectors, and take configuration sums over such states that are different from the ground state 
at finitely many edges. 

Let us restrict to the region h > 0. We conjecture that in the limit M, N — > oo the partition 
function ZjJ N of this lattice behaves as 

logZj^ ~ MiV (log /i«(C/e) + log /i (i) (C6) +iv(log^ ) (C) + log^ ) (-g- 1 r 1 ))- 

Here /z® is the partition function per site in the bulk theory, and i>W is the partition function per 
boundary site, which in the present normalization are given by 

//W(C) = 1 fori = 0,1, (3.1) 
^°)(C) = 1, ^)(C)=A(C;r). (3.2) 

These conjectures have been suggested by an argument similar to the inversion trick in the bulk 
theory, on the basis of the properties ( |2.8| )-( 2~?L3| ) and ( p. 22 ) for R and K matrices. Here we will 



not discuss the details. 

Now consider the lattice shown in Figure 3. We divide this lattice into the four sections indicated 
by the dotted lines. We denote the NW and SW corner transfer matrices (CTM's), by A^ w ((,^) 
and j4g|y(C£), and the upper and lower lines of K(C) matrices by \B)® and ^(B\. In the infinite 
lattice limit we can identify A^ w ((,^), Ag W ((,,£), \B)^ and as maps between, or elements 

of, the vectors spaces and (and their duals W*M and These vector spaces are 
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defined as the spans of the half-infinite pure tensor vectors 

■ • • ® «p(3) ® w p(2) ® u p (i), with = (-) j+ \ j > 1, for 

The identification is that 

Af w (C,0 ■■ H^^m, 

The partition function of the lattice is then given by 
3.2 Vertex Operators 

We define 'vertex operators' ^ 1_i,<) (0, 0* (1 ~ l>) (£), (CO and 0? ' (1_i,i) (C, by the half- 

infinite lattice insertions shown in Figure 4. The superscripts (1 — indicate how the vertex 
operators act on WW and WW. Namely, 

# (1 ~ M) (C,£) : n^^n^ </>? ,(1 ~ M) (c,e) m^n^. ( ' } 

Henceforth these superscripts will be suppressed. 

We now argue that the vertex operators ( |3.3|) and CTM's obey a set of relations that will allow 
us to use them to construct lattice realizations of 3> e (C)j ^e(C)> an d s(0| obeying relations 
(fO|), (§1| and (0). These relations are 

= 0I e (-g6> (3-4) 
E ^^(6/6)^(6)^(6) = <M6)<M£i), (3-5) 



e' e' 



E^(0^(C,0|fi) w = f w (O^(C,i/£0|£>« (3-6) 

e' 

£ w Wi?(c,i/0*$(0 = ^'WWfCCO, (3-7) 

e' 

4W(c,e)^(c,eo = &(e7o^(c,o, (3.8) 

41(c,e)C(evo = ^cco^cco, (3-9) 



where z/ l )(£') are given in (|3.2|). Properties ( |3.4[ ) and ( |3.5[ ) are consequences of the crossing symme- 
try of the i?-matrix and the Yang-Baxter equation (in the infinite lattice limit) respectively. They 



are discussed in reference [17]. The point in the argument for (|3.5| ) is as follows. If we compare the 



graphical expressions for both sides of (ft.5|), the left hand side contains one more site (or vertex), 
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Figure 4: The graphical definition of vertex operators 
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Figure 5: The boundary reflection property 
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Figure 6: The homogeneous CTM A^fa/fr) 



i.e., Relief (C1/C2)) m the formula, than the right hand site. The effect of the additional site is 
to multiply by the partition function per site /j". Since we normalized R so that [p® = 1, we 
have ( |3.5| ). Equation ( |3.6| ) follows from a similar graphical argument. See Figure 5. We use the 
boundary Yang Baxter equation in order to move the K(£') matrix up the boundary to infinity. 
This time we get the factor i/W . Equation ( |3.7| ) follows by a similar argument applied to the lower 
boundary. The argument leading to ( |3.8|) requires the introduction of the CTM of the homogeneous 
lattice denoted by A$ w (€) with one less argument. This is the CTM of the same vertex model 
defined on the homogeneous lattice of Figure 6. In [17| a simple graphical argument is given that 
leads to the equality 

^(0^(0 = ^(1)^(0- (3-io) 

This argument extends without change to the inhomogeneous lattice discussed here, giving 



4V(c,o<^(c,o = &(i)4k?(c,o- 



(3.11) 



Baxter has shown that the asymptotic (i.e. infinite lattice) limit of the CTM of the homogeneous 
8-vertex model is given by -<4.jvw(0 = /(0£°^> wnere is the CTM Hamiltonian (which is 



independent of an d /(0 is a scalar function. It is then a consequence of fljUOD that 

A%{t)UO = Me/0A { ^\0- (3-12) 
Finally, given (|3Tl2|) and property Q follows if 

4W(04V(C,0 = 9(C,t,?)A$ w ((,tf), (3.13) 



where 5(C)?)?') is some scalar function. In Appendix A we present a derivation of equation ( 3.13 ) 
using a generalization of Baxter's argument [16 for the asymptotic behavior of the CTM of the 
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homogeneous 8-vertex model. If we apply A$ w ((,£) to both side of (3.6), and then use (P-gp, we 



Similarly, 



obtain 

^^'(e')^(e7e)4W(c,oi5) w = ^(O^^^^oi^)®. 

Specializing to the case £ = 1, and making the identifications 

* e (C) ~ <MC), 

|0)fl ~ A^(C,l)|B)W and 
b<0| ~ (°)(S|4t(C,l), 

we obtain equations ( |2.3] ) and fl2.4| ) of Section 2. The remaining property ( |2.1| ) comes from ( p.5p . 

Thus, the vertex operators </> £ (£) and </>*(£), and the boundary states A^ w ((, 1)\B)^ and 
(°)(£|4?w(C> 1) are lattice realizations of the corresponding objects discussed in section 2. We 
speculate that the states \B)^ and ^{B\ will correspond to initial and final times states in the 



sense of [ 10 1 but now in an axial quantization scheme. 



3.3 Correlation functions and difference equations 



As in JTfl] , correlation functions of local operators of the inhomogeneous vertex model are defined 
in terms of 

(0 <B[^(c,o---^(c,o^^(c,e)A^(c,oC(c.o---^(c,e)iB>«, (3.15) 



where the superscripts on the CTM's are understood as modulo 2. From properties ( |3.8| ) and 
( |3~S| ) we see that it is possible to move all the vertex operators ^(£,£) and 0f (CjC) inside of the 
4vi^(C> £)^ivw(£' product. We can thus consider the more general spectral parameter dependent 
expression 

G (i HC,eieir--^) £l ,.,^ = (i) (5|4i(C,6^ 1 (6)-->^(e7v)4V(C0!^) (i) , with N even. 

(3.16) 

From Q3.14D G^°\(, £ = l|£i, • • • , £,n)ei,---,£ N is the lattice realization of the corresponding correlation 
function G(£i, ■ • • , £jv) E1 ,..., £JV of equation (|2~6|) . G (0) (C,£ = ■ • • ,£jv) E i,-, ejv thus obeys the 



difference equation (2.7) from the arguments already given in section 2. Here we shall attempt to 
give some insight into the physical origin of this difference equation by pointing out the sequence 
of steps necessary in order to derive it within the context of the lattice theory. For simplicity, we 
drop the e subscripts. At each step we refer to the vertex operator that arises as a result of the 
previous step. 



1 Move 0(£i) to the right of </>(£j+i) • • • 4>(£,n) using the exchange relation (|3^ 

2 Move </>(£-,-) to the right of the CTM Anw^, 1) using the commutation property Q3.g|). 
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kVl 

kVt 



3 'Reflect' </> u (C,£j) off the upper boundary using ( |376| 



4 Move <ft u (C, l/ij) back to the left of the CTM A NW ((, 1) using (j3j). 

5 Move <t>{l/Zj) to the left of 0(6) ■ ■ ■ 0(6-i)<K6+i) • • • <j>(£ N ) using Q. 

6 'Cross' (f)(l/£j) using the crossing relation (|3.4| ). 

7 Move 4>*{—q/£,j) to the left of ^4syi/(C> 1) using the commutation property (| 

8 Reflect 4> D ((, —q/£j) off the lower boundary using (|3.7|) , 

9 Move 4> D (C, -€j/q) to the right of A SW (C, 1) using (K% 
10 Cross <p*(—S,j/q) using (| 



11 Move 4>(Cj/q 2 ) to the right of 0(£i) • • • 0(£>-i) using (KB). 



This sequence of steps leads to the difference equation fl2.7|). One can of course derive the same 
equation by moving 4>{q~ 2 £,j) to the left of </>(£i) ' ' ' through 4sw(Cj 1) 3 reflecting off the 

lower boundary, etc. 

4 Boundary spontaneous magnetization 

In this section, we investigate the difference equation obeyed by the two point function of the XYZ 
spin chain. The solution of this equation, with a certain specialization of the spectral parameters, 
gives the boundary spontaneous magnetization. For general values of the boundary magnetic field, 
the solution does not factorize, but in the case of free boundary conditions, we find a factorized 
solution. As we show below, the value of the boundary spontaneous magnetization is (minus) the 
square of the spontaneous staggered polarization of the usual eight-vertex model without boundary. 

4.1 Two point functions 

Setting n = 2 in the elliptic version of ( 2.14T ), we have 

G(rti,C 2 ) = K 1 (Ci)« 2 i(CiC2)^i(Ci) J Ri2(Ci/C 2 )G(Ci,C2), (4.1) 

G((uq-%) = J R 2 l(^ 2 C2/Cl)^2(C2)i?12(ClC2)^2(C 2 )G(Cl, C 2 ), (4.2) 

where 

G(Ci,C 2 )= £ G *i.«(Ci>C2)v ei ®t> ea , 

ei,ea 

cWCi,C 2 ) = b(*1$ £1 (Ci)<MC2)K}b- 

In the following paragraph, we restrict our attention to the free boundary condition, r = — 1, 
or h = 0. The reflection matrices K and K in this limit are 



Here I denotes the 2x2 unit matrix and /(C) = f(C) ~ !)• 
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Let 

G±(Ci, C2) = G+-(Ci, C2) ± G_+(Ci, C2). (4.3) 
One can reduce equations (4.1), ( |4.2| ) to the following scalar difference equations: 



ft(rt,6) = ^#^ift(6,«. (4.4) 

/(Cl)/(-9 Ci) 

r /a _-v \ _ ^(£g" 2 C2/Ci)^(eCiC2) r , , , s 

W(Ci,<? C2) 777 \ 77 ZTT"^ — ^e(Cl,C2j, (4.5) 



where 

_ _i (-?C" 1 ;p)oo(-m" 1 C;?')c 



K (C 2 ) (-q(;p)oc{-pq X C ^pW 



and k(z) is defined in ( 2.18| ). We solve these equations by setting 

G e (Ci,C 2 ) = A(( 1 )A(C2)B £ (C 1 C2)B £ (Ci/C2)- 

The problem is now reduced to the following equations: 

A{q-\) 1 



A(C) /(C)/(-rt)' 

The solutions which are analytic in the region —q < £ < q~ 2 are given by 

(q 2 z~ 1 ;p, g 4 )oo(g 6 ^ _2 ;P, q 4 , q 8 )oo(pq 2 z~ 2 ;p, q 4 , q 8 ) 



A(C) 



B + (C) 



oc 



(pz 1 ;p,q 4 )oo(q 4 z 2 ;p,q 4 ,q 8 )oo(pq 4 z 2 ;p,q 4 ,q 8 )oo 
(q 4 z;p, q 4 )oo(q 10 z 2 ;p, q 4 , g 8 )oo(pq 6 z 2 ;iJ, q 4 , q 8 )^ 
(pq 2 z;p, q 4 ) (X (q 8 z 2 ;p, q 4 , q 8 )oo{pq 8 z 2 ;p, q 4 , q 8 )oo ' 
(-pq' 1 C 1 ;P,q 2 )oo(q 2 z' 1 ;p, q 4 ,q 4 )oo{pq 2 z~ 1 ;p,q 4 ,q 4 ) 



00 



(~q( 1 ;P,q 2 )oo(q 4 z 1 ;p,q 4 ,q 4 ) 00 (pz 1 ;p,q 4 ,q 4 ) o 
(-pq(;p, q 2 )oo(q 6 z;p, q 4 , q 4 )oo{pq & z]p, q 4 ,q 4 )oo 



{-q 3 (;p, q 2 )oo(q 8 z;p, q 4 ,q 4 )oo(pq 4 z;p, q 4 , q 4 ) 
B_(C) = 5+(-C), 



where z = 



4.2 Spontaneous magnetization 

The boundary magnetization is the vacuum expectation value of the boundary spin operator af . 
As in the XXZ chain, it is given by the value at £ = 1 of the following function: 

M (0) (Or) - B(0\E^(C,O-E W (t,Q\0) B 
U ' ' b(0\0)b 

where 

^5(Ci ! c 2 ) = ^ (i ' 1 - i) (Ci)^ M) (C2). 
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In terms of the functions (?e(Cl> C2) of (H? 



Using the result in the last section for the case r = — 1, we obtain 



M {0) (c -1) = (~ pz > p > q2 }° D (~ pz 1 ^P^ 2 )oo(-p;p,q 2 )lo 
(pz ; p, q 2 )oo (pz- 1 ; p, q 2 )oo (p; p, q 2 ^ 2 

^2„.„ „2\ ( n 2„-\. m „2\ f„2.„ „2\2 

I oc 



Joo 

(q 2 z;p, q 2 )oo(q 2 z~ 1 ;p, q 2 )oo(q 2 ;P, q 2 ) 



{-q 2 z;p,q 2 )oo(-q 2 z 1 \p,q 2 )oo(-q 2 ;p,q 2 )lo' 

The value at C, = 1 gives the spontaneous magnetization of the boundary spin operator: 

,2.„2\4 /_„. „\4 



w (o) n _ n _ (g ;g )5o(-p;p)5o 
l ' J " (-gW)^;^' 



(4.7) 

The XXZ limit p — > coincides with the previous result of |l3j . We again note the remarkable fact 
that the spontaneous magnetization -7W (0) (1; -1) is exactly the square of the corresponding bulk 
quantity, i.e. the spontaneous staggered polarization in the eight-vertex model [1^,17|. The same 
phenomenon was observed in the case of the XXZ chain in [13]. 

One can check formula (4.7) by comparing it with the derivative with respect to the external 
magnetic field h of the energy difference Ae(r) = eW(r) - e(°)(r). From pl28|) , we obtain 

dAe(r) _ 2 (q 2 ;q 2 )t (r 2 q 2 ;q 4 ) 00 (r- 2 q 2 ;q 4 ) 00 (p^p^^pf^r^p-p) 2 ^ 



9h (q 2 ;q 4 )lo( r q 2 ;q 2 )lo( r 1 q 2 ;q 2 )'L (p;p)io( r2 p;p 2 )oo(r 2 p;p 2 )oo' 

This quantity is equal to the difference of the magnetizations A4^(l; r) — A4^°'(l; r), which agrees, 
at r = — 1, with the result ( |4.7| ). 

One can also verify the result for directly using perturbation theory in e = 

—q/(l + q 2 )- As in [jl3"| , we solve order by order the equation 

E ( \(P*k+x°k + 1) + ) - 2e(Q + TQ') ] |0) B = 0, 
where 

k>l k>l 

The c-number normalization term Cfc(e) = Ylj>i c k,j £:i is included in order to ensure that the 
eigenvalue is zero. Solving for |0}_b, we find 

M {0 Hl; -1) = -1 + 8e 2 + 8(5r 2 - l)e 4 + 0(e 6 ). 

In terms of q and p, this expansion becomes 

A4(°)(l;-1) = -l + 8q 2 -24q 4 -8p + 0(q e ,pq 2 ,p 2 ), 

which agrees with the expression ( [4.7| ). 
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5 Boundary S matrix 



We now come to a brief discussion about the .S'-matrices which describe scatterings of quasi-particles 
with themselves or with the boundary. These S'-matrices obey the same type of Yang-Baxter 
equations, unitarity and crossing relations, as do the R- and the i^-matrices used to construct the 
model, but they are not the same. In fact, for the XYZ chain in the bulk, one expects [^] that the 
two-particle S-matrix has the same form as R but with the elliptic parameter p being changed: 

S(t; Pt q) = -R(S;p*,q), p* = pq' 2 . (5.1) 



A similar phenomenom was observed in |L3j], where the the boundary S-matrix M(£) for the XXZ 
chain was shown to be of the same form as the i^-matrix, up to an overall scalar factor, wherein 
the parameter r is changed to r* = rq^ 1 . 

We are then led to speculate that the M-matrix for the XYZ chain with a boundary (see (|6.2| ) 
for the characterization) is proportional to the iT-matrix, with both p and r scaled in the same 
way as above; namely 

M(£;p,q,r) = =— -K(£; P *,r*), p* = pq~ 2 , r* = rq' 1 , (5.2) 

f{^,p,q,r) 



where by K(£;p* ,r*) we mean the matrix K(£;p,r) = iT(£;r) in Q2.20) with the parameters p,r 



being replaced by p*,r* respectively. The scalar factor f(£;p, q, r) is fixed by solving the boundary 
unitarity and crossing relations, 

M(e)M(r x ) = i, 

M^-q-'r 1 ) = E s b : b - a a '(-qf)Mz;(o. 



with the condition that it reduces to the known result [13] in the XXZ limit p — > 0. Explicitly we 
have 



where 



f{£-,P,q,r) <p(£ 2 ]P*,q,r*)' 

lp(z;p*,q,r*) = <p (z; p* , q)?Pi(z; p* , q, r*) , 

* v (q 2 z 2 ;p*,q 8 )oo ( P *q 6 z 2 ;p*,q 8 )oz 

<Pi{z,P ,q,r ) 



(q 8 z 2 ;p*,q 8 )oo (p*z 2 ;p*, q 8 )^ 
(p*q 8 z 2 ; p* 2 , q 8 )oo(p* z 2 ; p* 2 , q 8 ) 



(p*q 4: z 2 ;p* 2 , q 8 ^ 2 



loo 



(gVz;/,? 4 )™ (p*r* 1 z;p*,q 4 ) o ( P *r*z;p* 2 ) 00 



(q 2 r*z;p*,q 4 ) OQ (p*q 2 r* l z;p* ,<? 4 )oc ( P *r* 1 z;p* 2 ) 00 



In the case of the XXZ chain, it was possible to derive the M-matrix since the boundary 
vacuum states are known explicitly in terms of the bosonic oscillators. Here we do not have such 



a construction, and (5.2) is no more than a plausible guess. We find it difficult to check it by 



18 



perturbative methods. As an alternative argument in its favor, we study below the continuum 
limit and show that the resulting formulas agree with those for the sine-Gordon model with a 
boundary [px}] . 

ice 

A 2A 



Introduce 



TT-2X 

and write 

f2 = (p^-Au/T, ^ = (p*f, r * = (p *)l/2-V7T_ (53) 

By the continuum limit we mean the limit K —* 00, K' — > tt/2, so that p* — > 1, A — > A, while 
keeping u and fj fixed. We find the following: 

limS(£) = p(u) x (sin(Avr) -(1 + a y <g> fj 27 ) + sin(A(vr - u)) ^(cj x ' (8) cr* + a z ® a z ) 

+ sin(Au) -(a x ® a x - a z ^cr 2 )), (5.4) 

IimM(0 = -Ro(u)a(0,u)a(fj,u)cos(Xu) x ( 1 - sm ( A ^ ] , (5.5) 

V cos 77 J 

Here the functions /o(u), Rq(u), <t(x,u) are given in (5.7), (5.21), (5.23) in [|i~0[| , respectively, and 
a(x,u) = a(x,u)/a(x,0), wherein our A and fj are to be identified with A and r\ there. f\ 



The formula (5.4) differs from the standard formula (5.6-7) [^] f° r the two-particle 5-matrix 
of the sine-Gordon theory by a gauge transformation. To see this, let 

i/= (J 7), tww, 

which has the property Uo x XJ~ x = o y , XJo y XJ~ x = a z , Ua z U~ l = a x . Let further 

V 

where U^' = U for j odd and = U for j even. Define the new basis of eigenstates by 

K». • • • >a)^,., m;W B = ^z£°te») • ••< 2) fe)< 1) (6)K) B - (5.6) 

In this basis the bulk and the boundary S'-matrices are given by 

S'(0 = (U®U) S(0 (u®uy X , M'(£) =UM(£)U~ 1 . 

Then lim S"(£) coincides with the formula (5.6-7) of [10[| . The limit of the M'-matrix is to be 
compared with a special case of the boundary S-matrix in [1C] 

| = 0, d = 0, 

where £ denotes the parameter £ in With this specialization the lim M'(£) agrees with (5.12), 
(5.21-25) in [ 10 1 up to an overall sign. In particular, the free boundary condition (5.29) [ 10 1 of 
the sine-Gordon theory is given by f] = tt(X + l)/2. It corresponds exactly to the free boundary 
condition r = — 1 in the XYZ chain. This is an indication that our speculative choice of the 
parameter r in the M-matrix (|5.2|) is exact. 



2 It seems that in (5.22), flOl, one should use a(x,u) = a(x,u)/a(x, 0) in place of a(x, u), to ensure Ri (0) = 1/ cos £. 
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6 Properties of boundary form factors 



We now consider the boundary form factors, defined as 

4° (&,..., 60= £ B(i\oK 1 &)---K n (tn)\i)B , (6.1) 



where O is some local operator, which commutes with type II vertex operators. Here, v* denote 
basis vectors in the space dual to V. The boundary form factor satisfies properties analogous to 
Smirnov's axioms for form factors, but with some differences due to the presence of a boundary. 
These properties follow from the properties of type II vertex operators and the boundary vacuum 
states. The commutation relations for type II vertex operators for the XYZ model were presented 
in [15]. In particular, 

KMi)K^2) = E ^(6)^(6)<\fkW6) , 



where the S-matrix is defined in ( |5.1[ ). 

By definition of the boundary S-matrix, type II vertex operators have the following property 
with respect to the boundary bound states: 

*:(0i<>fl = E JtfW «(0*j(r 1 )i»->B, 

e' 

same 1 ) = £B<i[*M<r 2 0^ i) s'(0. (6-2) 

where the boundary 5-matrices M (i) (£) = <r e AfW(-g -1 0<r B , M<W(£) = M(£) of Eq. Q, and 
(£; r) = cj^mW (£; r" 1 )^. 
Let 

(■■■®«* ®< »•••) ^-(o = £(---®<> 45(0 , 



e' e' 



5 < +1 (8 < ® • • ■) Pi, i+ i = ■ ■ ■ ® < ® < +1 ® • • • , 

(•••®<.®-.-)Mf (o = e(---®<'®---)m«:?(o, 



etc.. Then the boundary form factors defined in (|6.1| ) have the following properties, analogous to 
Smirnov's form factor axioms. 

Axiom I: Due to the exchange relation between type II vertex operators, which remains unchanged 
in the boundary theory, this property is identical to the usual case in the absence of a boundary: 

F n\Cl, • • • , £j, Cj+l) • • • £n) = -*n (£l> ■ • • > £j+l>£?> ■ ■ ■ £n) P ij'+l%j+l(£j/£?+l) . 

Axiom II: The analog of the periodicity condition is the difference equation satisfied by the form 
factors, which is similar to Cherednik's equation Q2.7|). It is a consequence of Axiom I and the 
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properties fl6,2| ) 

F®fa, q%, ...,£„) = • • • , • • • , ZuiSu-ifa/Sj-l) ■ ■ ■ SjA^/CiWfiq 2 ^) 

x S^tftitj) • • • 5 J --ij(g 2 ^--i^)S/ + i J -(9 2 6-+iCj) • • • SnjtfUj^tftj) 

Axiom III: The form factor ( |6.1| ) has simple annihilation poles due to the relation (see equation 
(B.9) in @0) 

* i I 1+61 



where "• • •" refers to regular terms and the scalar g* is as in Eq. ( p. 21 ) with p replaced by p* = p/q 2 ■ 
The residue of the form factor at the point £j = — g _1 £ri is 

ResF«(6, ■ ■ ■ • • • , = f i&(6, • • • , 6-1,0+1, • • • , Cn-l) ® C 

X Pn-2,n-l ■ • • -Fj'J+1 f<Sj,n-l(— £n/<Z£n-l) ' ' ' %j+l( — W^Cj'+l) 

- M n ((; n )Sj :n -i(-l/q(; n -it n ) ■ ■ ■ Sjj+ii-l/qZj+iZrJSjj-ii-l/qtj-iZn) 

x • • • S j , 1 (-l/qte n )Mf\-qZ- 1 )S 1>j (-qt 1 /Z n ) ■ ■ ■ S, qC, L /£ n )^ , (6.3) 

where C = £2 u e ® u -e- There are additional poles corresponding to bound states, which come 
from the poles in the S"-matrix and boundary S'-matrix. 

7 Discussion 

Let us summarize our results. We have shown that correlation functions of the semi- infinite XXZ 
spin chain obey the quantum Knizhnik-Zamolodchikov equation with reflection. This is a simple 
consequence of the exchange algebra of the vertex operators and the reflection property (|2.3j),( [2~4| ) 
of the boundary vacuum state and its dual. 

In the case of the XYZ chain, we assumed such a boundary vacuum state exists with a similar 



reflection property. Again using the exchange algebra between the elliptic vertex operators of [15] 



and the reflection property with the elliptic -matrix of [ 1 1 , 12 1 we were able to derive the q- 



difference equations satisfied by the correlation functions of the XYZ spin chain with a boundary. 
We cannot obtain an explicit formula for the vacuum states \i)g in this case, however, due to a 
lack of a bosonization formula for the elliptic vertex operators. 

The assumptions made in the latter case are supported by an alternative method of deriva- 
tion using the corner transfer matrix of the eight vertex model with a boundary corresponding 
to Sklyanin's transfer matrix. By constructing a lattice realization of type I vertex operators 
and boundary vacuum states, we showed that they obey the appropriate commutation relations 



3 There is an error in equation (B.9): The first factor is the inverse of the correct expression. 
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and reflection condition, and therefore that the correlation functions obey the correct difference 
equations. Using the difference equations, we computed the boundary spontaneous magnetization 
M(°)(l;-1), i.e., the magnetization in the case where the boundary field vanishes. 

We also conjectured a natural form for the boundary S-matrix for the XYZ model, which is the 
reflection matrix associated with type II vertex operators. In the continuum limit the boundary 
S-matrix becomes (up to a gauge transformation) that of the boundary sine-Gordon model JH| . 
Finally, using this boundary S'-matrix and the exchange algebra of the elliptic type II vertex 
operators, we formulated the analogues of Smirnov's axioms for form factors H. 

For the XYZ chain, neither the bulk nor the boundary S'-matrix have been derived using an 
alternative method to the vertex operator approach. It should be possible to obtain these using 
the Bethe ansatz approach. The same can be said for the boundary S-matrix of the XXZ chain. 

In our boundary S-matrix, one can find a pole £ 2 = rq^ 1 in the physical strip 1 < |£ 2 | < \q\~ 2 - 
This pole yields a boundary bound state, and as in the XXZ chain case |l3| this state should be 
identified with the second vacuum state \1)b- The coincidence of the vanishing point r = q, —q of 
dAe/dh and the point where the pole leaves the physical strip gives a consistency check for the 
conjectured form of the boundary S-matrix. 

In section 6, we derived the difference equation for form factors of the semi-infinite spin chain. 
In the continuum limit, one can derive the difference equations satisfied by form factors in massive 
integrable quantum field theory on a semi-infinite line. For example, for the sine-Gordon theory, 
using the setting given in Section 5 and taking the limit p* — > 1 of the equation ( |6.3D , one finds 

Fniftir ••,/% + 27TZ, • • • ,/3 n ) 

xSij((3j + Pi - 2-ki) ■ • • Sj-ij(j3j + P j+ i - 2m) 

xSj+ijiPj + Pj-i ~ M • ■ ■ S nj (pj +Pi- 2m)Mj(Pj - 2m) 

xS jn (pj -p n - 2iri) ■ ■ ■ S jj+ i(Pj - p j+ i - 2m), (7.1) 

where we also set u = —iP and 

]imF^i,---,^ n )=^ n {Pi,---,P n ), 

limS'(0 = S(P), limM'(e) = M{fi). (7.2) 

Here the primed quantities are the gauge transformed ones (see Sec. 5). To solve this equation 
is an open problem. It seems possible by applying the method developed by Smirnov [|| with 
modification by Sklyanin's Bethe ansatz scheme ||. 

Taking the quasi-classical limit of our difference equation, one can derive the boundary analog 
of the KZ equation. It is however unclear whether this equation governs the boundary conformal 
field theory in the way that the KZ equation does for the bulk theory. It would be nice if one could 
relate the equation to the one obtained by Cardy Q. 
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A The Asymptotic Behaviour of the CTM 



In this appendix we derive the asymptotic form of the CTM of the inhomogeneous 8- vertex model 
using a generalisation of Baxter's argument for the homogeneous model. Consider the transfer 
matrix T({Q}, shown in Figure Al. 



Ci (2 ■ ■ ■ Cn-i (n 



6 


6 










£jv-i 

























Figure A.l: The finite completely inhomogeneous transfer matrix T({Q}, 

It is the trace of the product of N R-matrices, each with independent vertical and horizontal 
rapidities Q and i = 1,---,N. Two such transfer matrices T({Q},{^i}) and T ({£[)-, {^}) 
commute if (£i/Ci)/(£i/CD =const, independent of i. This means that the normalised eigenvectors 
of T({d}, {&}) depend only on the ratios (£i/Ci )/(£?/£?) for z ^ j @. 

The first step is to specialise the transfer matrix to the case when all the horizontal rapidities 
on the left half of the lattice are equal to the left vertical rapidities are alternately £' and 1/C 
the right horizontal rapidities are all £ and the right vertical rapidities are alternately £ and l/£. 
Then we identify the infinite product of such transfer matrices (taking N — > 00 as well) with the 
product Asw(Ci €')Anw(Ci °f suitably normalised CTM's. Such an identification implies 

AswiC, a')A NW ((, = r(c', c, e, o*(C, c, e/o- ( A -i) 

Here and elsewhere we adopt the convention that lower case letters represent scalar functions and 
that upper case letters represent matrix functions. We also suppress the i superscripts on the 
CTM's. Now send £ — ► x/£' an d £' — ► x/£ (where x is an arbitrary constant) and eliminate X. 
Then, 

r(C', C, x/£, x/O^WC, O^W(C, = t(C', C, £)^W(C, X/0A NW (£, x/O- (A.2) 

This equation immediately tells us that the product 

Anw{(, = A NW (C, OAnw(C, (A.3) 

(where fi is another constant) depends on £ only through a scalar function. Setting t; = fi in QA.2| ), 
solving for AswCC'jO and substituting in (|A.1|) gives 



x/0^w(C = c(C, C, OnC, C, £/0> (A.4) 
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where 

nC,C,£/0 = ^(Cx/M^CCC.e/O^wCC/i). 

(We shall suppress dependencies on the constants x an d A 4 -) Setting C = C = 1, we can rewrite 
©as 

A NW (1, S')A NW (1, = dtf, t)Zg£), (A.6) 

where = c(l, 1, £) and = F(l, 1, We sha11 solve this equation for 
AnwO-,0 an d then reclaim 

A NW ((,0=r((,OA N w(l,0, (A.7) 
where r(£,£) is the scalar function r(£, £) = Anw(^,^)A]^ w (1,^). 



Equation ( |A.6[ ) is the generalisation of (13.5.17) of Baxter's book |)lq| , and the argument now 
proceeds as in [O] with some minor modifications. Interchanging f and £ and eliminating Z gives, 



Now consider a representation in which ^4jvw(l)C) is diagonal. Then ( |A.8| ) implies 

= , A9) 

Thus Ajv^(1,^'), .Ajvw(1>£) and Z(£'£) commute and have common eigenvectors, independent of 
£' and £. We can clearly renormalise in the following way: 



A* w (i,0 = A NW (i,0/a(0 
z R {t) = z(0/P(0, 



„,.w„,„ (A. 10) 



where a(£) and /?(£) are the respective eigenvalues of some common eigenvector of AnwO-iQ an d 
Z(£). Then equation (|A.6|) becomes 



AK w (Z')A% w (0 = Z R m- (A.11) 

Now we can diagonalise, 

A d NW (0 = P- 1 A* w (OP, (A.12) 
where P is the matrix of eigenvectors of A^ w (^). Define Z d '(£'£) similarly such that 

A d NW (OA d NW (0 = Z d m- (A.13) 
Differentiating with respect to £, it is apparent that the general solution is 

A d NWrr (0 = m r C ar r = 1, 2, • • • , (A.14) 



where m r and a r are independent of £ and £. Using the redefinitions (A. 3), ( A.7|) , ( |A.1C| ) and 
( |Al2D , we find, 



Anw(C,0 = a(C,OPA d NW (OQ~HO, (A-15) 
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where Q = P 1 -Ajvw(C) A 4 ) and the scalar function a(C£) = r (C>£) a (£)- Setting £ = 1 in 
( [A. 15 ), solving for Q _1 (C) an d substituting back into the same equation gives 

A NW ((,0 = ^^ ) PA d NW (OA d NW (l)- 1 p- 1 A NW (C,l). (A.16) 

Defining the operator D d as that with diagonal entries a±,a2, ■ ■ •, gives 

Anw((,0 = ^\r D A NW (C,l), (A.17) 

where £~ D = PC^P' 1 - 

We may now specialise this formula to the case when £ = 1 and all the horizontal arrows of 
Anw point to the left (the argument is independent of the direction of the arrows). The CTM is 
that of Figure 4, and is denoted by Anw{£,)- Since Anw{^) = id, this gives 

A NW (0 = ?^r D (A.18) 
o(l,l) 

where D is unchanged and the bar indicates that the scalar factors are those relevant to the 
homogeneous CTM Aj^wiO- Hence we obtain the desired result Q3.13 ) that 

A NW (OANw(C,e) = g((,Z,e)A NW ((,tt') where, 

a(C en - «(C,ffi(l,fl (A. 19) 
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